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Abstract. With future wide and deep cosmological sky surveys, a large number of 
gravitationally lensed, multiply imaged systems will be found. In addition to multiply 
imaged galaxies and quasars, sources will include transient events like supernovae and 
gamma ray bursts in which case very accurate time delay measurements are possible. 
Also, large numbers of systems with several lensed sources behind a single lens will 
be observed. In this paper, we review and compare different possibilities of using 
future strong lensing data to probe lens matter distributions and to determine the 
Hubble parameter and the matter density of the universe. Specifically, we investigate 
the possibility to break the well-known degeneracy between dark matter halo profiles 
and the Hubble parameter using observed flux ratios. We also investigate how 
strong lensing can provide useful constraints on the matter density of the universe 
independently of the flux ratios and other cosmological probes. 



1. Introduction 

In the last decade, a picture of our universe being dominated by dark energy and dark 
matter has emerged. Elements of this picture include observations of Type la supernova 
(SNIa) distances [1-4], anisotropies in the cosmic microwave background (CMB) [5-10] 
and the large scale structure (LSS) of galaxies [11, 12]. For a recent concise review of 
cosmological bounds on dark matter and dark energy, see Ref. [13]. A flat universe with 

~ 0.7 and ~ 0.3 provides a good fit to all current observations, indicating a 
recent {z ~ 0.5) transition from a decelerating universal expansion to an accelerating 
universal expansion. The current expansion rate (and hence overall scale) of the universe 
as measured by the Hubble Space Telescope (HST) Key Project [14] is h = 0.72 ± 0.08 
where h is the dimensionless Hubble constant, h = /fo/(100kms~^ Mpc~^). 

Though the picture described above is very effective in explaining all observational 
data with just a few free parameters, the nature of the dark energy and dark matter is 
still unknown. The simplest dark energy model is the cosmological constant. Difficulties 
in theoretically explaining the magnitude of the energy density has brought about a large 
number of alternative explanations, including quintessence [15-17], k-essence [18, 19] and 
phantom energy [20,21]. However, current observations are still best explained by the 
cosmological constant [22,23]. 
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On large scales, dark matter is well- described as being non-relativistic and 
collisionless, i.e. cold dark matter (CDM). On galaxy scales however, there has been 
a dispute whether there is a discrepancy between numerical simulations of galaxy 
formation in CDM models [24-28] and observations of rotation curves of dark matter 
dominated galaxies [29-31]. Precise observational determinations of halo density profiles 
will help in understanding whether there is any contradiction with the simulations and 
provide further insight into the properties of dark matter. 

When investigating the matter distribution of the universe, gravitational lensing 
has the advantage of being equally sensitive to all kinds of matter, regardless of it's 
microscopic properties (e.g. being baryonic or non-baryonic) and dynamical state. Both 
strong lensing (multiply imaged sources) and weak lensing (weakly distorted source 
images) has been used to constrain the matter distribution in galaxy scale [32-34] and 
cluster scale [35-38] halos. 

Besides being able to provide information on the matter distribution of dark matter 
halos, gravitational lensing has also the potential to constrain cosmological parameters 
such as JIm, and h. More than 40 years ago, Refsdal showed how the time delay 
between multiple images of supernovae (SNe) could be used to measure h and the mass 
of the lensing galaxy [39] . To date - for galaxy mass lenses - approximately 100 multiply 
imaged sources (including galaxies and quasars) have been observed out of which there 
are eleven well-determined time delay measurements (all quasars) [32, 40, 41]. For several 
reasons, it has been difficult to implement Rcfsdal's method. First, time delays between 
quasar images arc notoriously difficult to measure compared to what would be the case 
for transient sources like supernovae. Second, there is a degeneracy between the derived 
value of h and the matter distribution of the lensing galaxy. Modeling the lenses as 
more concentrated gives larger values of h and vice versa [42, 43]. 

Weak lensing can be used to probe the matter distribution of the universe also on 
very large scales. For reviews on the current state of cosmological constraints from LSS 
weak lensing, see Refs. [44-46]. There has also been attempts to constrain Vt^ through 
the statistics of gravitationally lensed quasars with disparate results [47-50]. 

Future strong lensing data will include observations from e.g. Pan-STARRS|, the 
Large Synoptic Survey Telescope (LSST)§, the Joint Dark Energy Mission (JDEM) 
contender the Supernova Acceleration Probe (SNAP)|| and the James Webb Space 
Telescope (JWST)^. In this note, we use simulated strong lensing data for a SNAP 
like mission combined with follow up observations from, e.g. JWST to investigate 
and compare different possibilities of constraining the slope of lens density profiles, the 
Hubble parameter and the matter density of the universe. 

In Sec. 2 we present our lens model. In Sec. 3, we investigate different uses of a data 
set consisting of a large number of lens systems, each with a multiply imaged source 

\ http://www.pan-starrs.org 

§ http://www.lsst.org 

II http://snap.lbl.gov 

^ http : //www . j wst . nasa . gov 
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with well-determined time delay. This section extends and generalizes earlier work in 
Refs. [51] and [52] (sec also Refs. [53-56]). In Sec. 4, we compare our results with what 
can be obtained by studying lens systems with several multiply lensed sources in a single 
lens (cf. the case of Abell 1689 [36]), with or without measured time delays. Finally, 
Sec. 5 contains a summary and a discussion of our results. 

The cosmology used in the simulations throughout this paper \sh — 0.65, Vtyi = 0.3 
and JIa = 0.7. Our default dark matter halo is the singular isothermal sphere described 
in Sec. 2. 



2. Lens model 

The lens equation relating the angular source position 9 and the angular image position 
/? is given by [57] 

pDs = ODs - dDis, (1) 

where a is the light deflection angle and Dig, Di and Dg arc angular diameter distances 
between source and lens, lens and observer, and source and observer, respectively. The 
lens equation can be rewritten in dimcnsionlcss form by introducing an arbitrary length 
scale ,^0 in the lens plane and a corresponding length scale rjo = ^oDg/Dd in the source 
plane. In terms of the dimensionless vectors x = 9Di/^q and y = PDs/rjo the lens 
equation takes the form 

y = x-a{x), (2) 

where 

a{x) = ^^^a{iox) (3) 

is the scaled dcfelction angle. 

In this paper, we assume that the matter distribution in lenses can be described by 
a simple power-law density profile 

p(r) oc r~'^. (4) 

Since almost all images of multiply lensed sources are at a limited range of small r, 
this is a good approximation despite the fact that the slope is expected to change with 
radius at large r [24-27]. Note that the baryonic contribution to the density profile of 
lensing galaxies is non-negligible at small radii and that results from strong lensing not 
necessarily reflect the nature of pure dark matter halos. 

For 77 = 2, we obtain the familiar singular isothermal sphere (SIS) model built upon 
the assumption that the mass components behave like particles in an ideal gas, confined 
by their spherically symmetric gravitational potential; 

Here, v is the line-of-sight velocity dispersion of the mass particles. 
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For all lensing purposes in this paper, matter distributions can be described by the 
projected Newtonian gravitational potential 



(6) 



For the power-law density profile, the projected potential is given by [40] 

* (X x^-" (7) 

where x is the impact parameter (in arbitrary units, ^o), i-e. the minimum distance 
between the light ray and the lens centre. 



2.1. The Einstein radius 



Since the scaled deflection angle a = V^, for a spherically symmetric lens we have 

* = ^/-'^ (8) 

where a;^ is the Einstein radius; the solution to the lens equation for y = Q. The lens 
equation for the power-law lens can then be written (where x > and a is directed 
towards the center of the lens) 

(9) 



y — X — a — X — x^ 
For multiple images, we have 

y — xi — ai — a2 — X2 



and 
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Putting = Di, i.e. denoting positions in terms of angles, we get 



Of 



'7-1 



01 



SIS, 



where 6'sis is the Einstein radius for a SIS halo, 

'v\^D, 



(10) 
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2.2. Time delay 

The time delay for a gravitationally lensed image as compared to an undeflected image 
is in the general case given by [57] 



T 



D,D, 



;i+^o 



(x-y) 
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(14) 



The time delay between two lensed images is given by At — T2 — Ti. For an isothermal 
lens, we have 



At 



SIS — 7: 



2 D, 



Is 



[l + zi){ei-e',)^2aDi{l + zi)p. 



(15) 
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Note that the time delay is independent of the source redshift. In the general case of 
T) ^ 2, we can rewrite Eq. (14) in terms of g = A6/ <6 > where A^^ = 9i — 62 and 
<9>— {9i + ^2)72 and Taylor expand to obtain (see also Refs. [40] and [58]) 

{r] - 1) DiDs 



At 



D, 



~{r)- l)At 



SIS 



1 - 



(2 - Vf 
12 



(16) 



There is a strong degeneracy between rj—l and h (which comes in through the distances). 
This degeneracy can be broken by including data from the observed image luminosities; 
the flux ratio. 



2.3. Flux ratio 



We denote the flux ratio r = |/ii//i.2|. The magnification /i for a source at position y 
observed at position x is given by 

X dx 



y dy 



(17) 



For a power-law lens, we can express the flux ratio as 



= rsis 



[2-r^)9l-\9, + 92)/{9l-^ + 9l-^) 



l-{2-7^)9l-\9, + 92)/{9i-'' + 9 



2-ri 



l^(2-r^)9l-\9, + 92)/i9t-'' + 9 
Expanding this to second order in q, we obtain 

r ~ 1 + (7^ - l)g + ^(7^ - ifq^ + 0{q^). 
Note that the flux ratio is independent of the value of h. 



(18) 



(19) 



3. Individual multiple image systems 

With a SNAP like JDEM contender, a large number of core collapse supernovae (CC 
SNe) will be discovered [59]. Monte Carlo simulations performed using the SNOC 
package [60] predicts that out of a total of ~ 10^ CC SNe at 2; < 5 in a 20 square degree 
field during three years [61], ~ 850 will be multiply imaged with an I-band or J-band 
peak brightness < 28.5 mag for the dimmest image+. We demand that the surface 
brightness of the lens galaxy is fainter than 24 I magnitudes per square arcsecond at 
the position of the faintest image, in order to avoid contamination from the lens galaxy. 
In accordance with Ref. [52], we also introduce a quality factor, /, giving the fraction 
of the lens systems that are "simple", i.e. a single dominant deflector with moderate 

+ Since these simulations were performed, the SNAP specifications has been shghtly degraded and the 
corresponding point-source magnitude Umit is now 27.7. This number would yield ~ 400 multiply 
imaged CC SNe. 
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ellipticity and external shear. In the following, we assume a quality factor, / = 0.5, 
in accordance with currently observed systems [40]. This cut leaves a total of ~ 400 
systems. 

Here, we investigate how measurements of the image positions, time delays, flux 
ratios and lens galaxy properties for such systems could be used to constrain the slope of 
lens density profiles, 77, the Hubble parameter, /i, and the matter density of the universe, 

3.1. Error budget 

Our error assumptions basically follow Refs. [59] and [52]. We conservatively use a 
positional uncertainty of ae — 0.01". Combining data points from several filters could 
be combined to give a time delay uncertainty (jAt of 0.05 days for the SNAP mission. 
However, this number may be degraded if the lightcurve is modified by microlensing. In 
the following we set dAt = 0.15 days. For the flux ratios, wc assume the error ar to be 
dominated by microlensing [62,63] and lensing by CDM substructure [64-66]. An error 
estimate of Ur/r = 0.5 is employed. Ideally, one would like to use the proper probability 
distribution functions for the microlensing magnification probability as a function of, 
e.g. position in lens galaxy. In the following however, we assume a simple gaussian 
distribution and assume that extreme events with very large microlensing magnification 
factors can be identified from, e.g. pecularities of the lightcurves and be removed from 
the sample. This will however only be true if the microlensing timescale is smaller 
than the timescale of the intrinsic source variations. Redshifts error are negligible; we 
adopt (72; = 0.001 from photometric redshift measurements based on SNAP multi-band 
photometry. 

Even though we can derive lensing statistics using simple spherical galaxy profiles, 
it is important to include lens ellipticities and galaxy environments in the analysis [67]. 
Motivated by models of known gravitational lens systems, we assume an uncertainty 
in the alignment between the mass distribution and the optical light distribution of 
(A0^)^/^ < 10 degrees [68]. Prom N-body simulations and semianalytic models of galaxy 
formation, we can estimate the levels of external shear due to structure near the lens 
in gravitational lens systems to be 7ext = 0.058, with an rms dispersion of 0.071 [69]. 
We modify o\ii simple simulated lensing systems to include this effect. We also add an 
uncertainty due to scatter in the slope of individual halos of Ar] = 0.2. 

The addition of a constant surface density k^. to a lens, changes the time delay by a 
factor 1 — Kc but leaves the image positions and flux ratios unchanged. In the following 
we have assumed the effects from external convergence to be comparable to or smaller 
than the effects from external shear and a scatter in r). However, it should be noted 
that the effect potentially can be very important [70, 71]. The error in the derived value 
of h will be directly proportional to the systematic error in the assumed value of the 
external shear. 

When constraining h, we use a prior on the matter density of 0,^ — 0.30 ± 0.04 
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Figure 1. Constraints in the [77, /i]-plane using ~ 400 multiply imaged SNe 
(corresponding to a quality factor / = 0.5). In the left panel, constraints from the 
time delays are shown, in the middle panel constraints from the flux ratios. Combined 
constraints are shown in the right panel. Contours correspond to 68.3%, 90%, 95% 
and 99 % confidence levels for two parameters. The black line indicates the 95 % 
confidence level (2cr) for one parameter. 



[11]. Results are very insensitve to the exact value and size of the error on this prior. 
3.2. The Hubble parameter 

For constraints in the [rj, /ij-plane, we follow the same recipe as Ref. [52] and obtain very 
similar results, see Fig. 1. Small differences are due to the random fluctuations in the 
simulation of lens systems as well as the fact that we only include the scatter in the slope 
T] in the uncertainty - not as a modiflcation of the simulated sample as was the case in 
Ref. [52]. This is to ensure that our confldence contours are centred at the correct value 
to facilitate estimates of the contour sizes. Constraints from the time delays are shown in 
the left panel and constraints from the flux ratios in the middle panel. In deriving these 
constraints, we have used Eqns. (16) and (19) generalized to include the effects from 
galaxy ellipticities and external shear, see Appendix A. Contours correspond to 68.3 %, 
90 %, 95 % and 99 % confldence levels for both parameters to have values within their 
borders. The black line indicates the 95 % confldence level for one of the parameters 
to lie within the contour (i.e. 2a). Because of the strong degeneracy in the time delay 
between h and r], flux ratio measurements are needed to constrain the Hubble parameter. 
Combining the results from the time delay and flux ratio measurements, we are able to 
make a determination of h within 10 % and - perhaps more interestingly - to determine r) 
at the per cent level at 95 % confldence (right panel). Even for an extremely conservative 
value of / = 0.05, we can determine t] within 2 % and h within 25 %. It should be noted 
however that these numbers may be severly degraded in the presence of systems with 
very large microlens magniflcations. 
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Figure 2. Constraints in the [77, r^Mj-plane using 400 multiply imaged SNe 
(corresponding to a quality factor / = 0.5). In the left panel, constraints from the 
time delay ratios are shown, in the middle panel constraints from the flux ratios. 
Combined constraints are shown in the right panel. 



3.3. The matter density 

Since the time delay is proportional to h through the distances, we can study the ratio 
of time delays in order to factor out the Hubble constant and constrain (assuming a 
flat universe). I doing this, we wish to study ratios of lens systems with large differences 
in redshift in order to maximize the sensitivity of the ratio to the value of the matter 
density. Denoting ki = At^/At^, where m and n numbers the lens systems in the total 
sample, we want to combine the time delays (i.e. choose m and n) in order to maximize 
{dk / dVLy[) / ak = k'/a^ (i.e. maximizing the sensitivity to flu)- This is equivalent to 
maximizing the quantity 

At'JAtra - At;/At„ 

y^AtJTAt^ 

From a sample of N lens systems, we are thus able to obtain — 1 time delay ratios 
with maximum sensitivity to Qm- Note the similarity of the method to the use of cross- 
correlation tomography in weak lensing [72]. The time delay is thus used to constrain 
ratios of the distance combination (DiDg)/ Di^, see Eq. (16). The flux ratio data is used 
to constrain rj as described in Sec. 3.2. In Fig. 2, results from the time delays (left panel), 
flux ratios (middle panel) and the combined results (right panel) are shown. Note that 
the time delay ratios basically constrain Qm independent of the value of rj which allows 
us to constrain the matter density independently of the flux ratio measurements that 
may be plagued by large uncertainties. The "false" minima at very low flu is due to 
the fact that the ratio of distances often is degenerate in Qu- The same effect can be 
seen in Fig. 3 that basically probes the distance combination Dig/Dg, see Eqns. (12) and 
(13). Though not comparable in precision to other probes of the matter density such 
as the LSS of galaxies, strong lensing should be able to rule out an Einstein-de Sitter 
universe with = 1 at very high confidence level and provide an independent sanity 
check of the matter density. Setting / = 0.05 leaves more or less unconstrained. 
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4. Multiple source systems 

We next compare our results obtained with single source systems with the possibilty to 
constrain the matter density in the universe and the matter distribution in lenses using 
systems with multiple sources being lensed by the same lens. The cross-section for 
multiple imaging is larger for cluster mass lenses than for galaxy size halos. However, 
the time delays for cluster lenses are generally very large. Also, if the sources are 
galaxies, we arc not able to measure the time delays. Examples of such systems are 
Abell 2218 [74] and Abell 1689 where 30 background galaxies are being lensed by the 
cluster lens and displays more than 100 images [36]. 

The cross-section for multiple imaging by a SIS halo is given by 



where r is the flux ratio. For a lens redshift of zi — 0.5 and a source redshift Zg — 1.5, 
we get 



^.1. No time delay measurements 

Since the cross-section for multiple imaging is proportional to f^, the probability for 
multiple sources to be strongly lensed by a galaxy size halo is small. However, in 
sufficiently deep exposures of any given massive clusters, we expect to find a number of 
sources being lensed into multiple images. The vast majority of these sources will be 
galaxies in which case we will not be able to measure the time delays between images. 
Also, for massive cluster size lenses, time delays are generally too long to be practically 
observable even for variable sources. [The cross-section for two images with brightness 
ratio less than r is oc (^) . Thus, we expect to observe relatively few systems with 
r ~ 1 and small time delays At, see Eq. (23).] 

The Einstein radius is a function of both the halo slope and cosmological distances, 
see Eqns. (12) and (13) and in cases where we have multiple images of sources at different 
redshifts, we can use the observed image positions to constrain rj and VLm- A sample of 
28 multiply lensed galaxies with redshift measurements or reliable photometric redshifts 
behind Abell 1689 gives a constraint VLy[ + < 1.2, i.e. the data is in accordance with 
an Einstein-de Sitter universe with VLm = 1- The constraining power of the sample is 
limited by the low redshift {z ~ 0.18) of the cluster [73]. However, an analysis of four 
multiple image systems in Abell 2218 (also at 2; ~ 0.18) gives < Qm < 0.3 assuming a 
flat universe, ruling out an Einstein-de Sitter universe at 5a [74] . The difference in the 
claimed accuracy of the results warrants further analysis of the sensitivity of the method 




(21) 



and the time delay for the two images is 




(22) 




(23) 
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Figure 3. Results from fitting the image positions for 16 sources at regular intervals 
between z = 1 and z = 4.75 behind a SIS cluster at z = 0.5 with v — 1000 km/s. The 
uncertainty in the oberved image positions is ae ~ 0.2". 



on the model assumptions and uncertainties. This question has been investigated in 
great detail in, e.g. Ref. [75]. Here, we investigate the sensitivity of the bound on Qu 
on the slope of the dark matter halo and the size of the observational uncertainties. We 
simulate a gravitational lens system consisting of a massive SIS cluster at z = 0.5 with 
V = 1000 km/s. We then distribute 16 sources at regular intervals between z = 1 and 
z = 4.75 with random source positions (inside the cross-section for multiple imaging). 
Note that this analysis is very primitive in the sense that we have not included any of 
the complexities of real cluster lenses such as substructure and departures from spherical 
symmetry. However, this simplicity allows us to compare the efficiency of the method 
compared to the ones described in Sec. 3. The accuracy of our results is completely 
determined by the observational error in image position. In Fig. 3, results for ag = 0.2" 
is shown, yielding error contours comparable to the ones in Ref. [74] for ^m- Setting 
ag = 0.5" basically leaves flu unconstrained while ag = 0.02" yields = 0.3 ± 0.04 
and constrains t] at the per mille level*. Results for rj are very sensitive to the exact 
source positions since large values of r dramatically improve the limits on rj, see Sec. 5. 

Our analysis shows that for a massive cluster with of the order 20 multiply imaged 
background sources at a large redshift range, we need the image positions determined 
to an accuracy of better than ag ~ 0.2" to give useful bounds on ^m- This number 
can be relaxed if the number of sources is larger. The result is insensitive to any prior 
information on rj and thus on the error in the observed flux ratios (assumed here to be 
a^/r = 0.5). 

4-2. Time delay measurements 

As explained in Sec. 4.1, we expect to observe very few lens systems with multiple images 
of several variable or transient sources with time delays short enough to be practically 
measureable. Nevertheless, we here explore possible constraints from such systems. In 

* Of course, at this level of precision, the lens model is too crude. 
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Figure 4. Results from the time delay ratio (left panel), flux ratios (middle panel) and 
the combined results (right panel) for a SIS lens at zi = 0.5 with a velocity dispersion 
oi V = 300 km/s and two multiple imaged quasars at Zs,a = 1 and Zs^b = 3 at angular 
source positions (3a = 0.5" and = 1.0". 

order to have a non-negligible cross-section for multiple imaging but reasonable time 
delays, we imagine lensing from a very heavy galaxy or a small group of galaxies. Our 
default system is a SIS lensing halo at zi = 0.5 with a velocity dispersion of v = 300 
km/s. We imagine two multiply imaged quasars at Zs,a = 1 and Zs^b = 3 at angular 
source positions Pa = 0.5" and Pb = 1-0". The assumed errors are the same as in 
Sec. 3.1 except that we set the error in the time delays to a At = 2.0 days, corresponding 
to typical errors for time delays between quasar images [40]. For such a configuration, 
we have a time delay for the images of source "a" of Ata = 148 days and a flux ratio of 
Ta = 2.6. The corresponding numbers for source "b" is Atb = 296 days and = 3.4. 
We derive constraints on rj and Qu from the observed flux ratios and the time delay 
ratio (c.f. Sec. 3.3), see Fig. 4. The error budget is equally shared between the error 
in the time delay measurements and the observational error in the image positions of 
ag = 0.01". Therefore does our results not improve significantly for transient sources like 
SNe or gamma ray bursts where more accurate time delay measurements are possible, 
unless (7g can be significantly reduced. Note that the shape and size of the contours in 
the [rj, f2M]-plane is very dependent on the exact configuration of the lens system. In 
general, results benefit from a large redshift separation for the sources since the time 
delay ratio is proportional to {Dig/ Ds)a/ {Dis/ Ds)b- However, even in the lucky case 
of finding a lensing system with two well-observed multiply imaged sources with time 
delays determined to good accuracy, we can only hope for very rough constraints on r) 
and flM- 

5. Summary and discussion 

We have analyzed and compared how measurements of the image positions, flux 
ratios and time delays in different types of strong lensing systems will allow for the 
determination of the matter distribution in galaxies, as well as global cosmological 
parameters such as the Hubble constant h and the matter density ^m- While CMB 
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and LSS data probe properties of dark matter on large scales, measurements of galactic 
halos can probe the small scale properties of dark matter. Furthermore, gravitational 
lensing mainly probes halos at relatively high redshift where few methods for measuring 
galaxy density profiles are available. 

First, we investigated different uses of a data set consisting of a large number of 
lens systems, each with a multiply imaged source with the time delay measured to high 
precision. Such a data set can constrain the slope of dark matter halos 77 at the per cent 
level and provides an independent test of the values of h and JIm- Combining the results 
from the time delay and flux ratio measurements, we are able to make a determination 
of h within 10% and 7/ within 1% {2a). 

Factoring out the Hubble constant by studying time delay ratios, we are able to rule 
out = 1 at high confidence (assuming a flat universe). This bound is independent 
of the slope of the halos and the flux ratio measurements. 

Second, we compared the results obtained for multiple systems with single sources 
with the use of lens systems with several multiply lensed sources in a single lens. For 
massive cluster size lenses, we expect to find several multiply imaged galaxy sources. 
The position of these images can provide a sanity check on Qm- The error budget is 
dominated by the observational error in image position that need to be of the order or 
less than ~ 0.2" to be able to rule out an Einstein-de Sitter universe. Systems with 
multiply imaged sources with measurable time delays will be rare and not very useful 
for the purposes discussed in this paper. 

In the case of a large sample of individual multiple image systems, the error budget 
in the time delay analysis is dominated by the observational time delay error. The flux 
ratio constraints are quite robust to an increase in the (gaussian) error size. Increasing 
the errors by a factor of 3 (to 150%), we are still able to obtain 77 = 2 ±0.05 (2(t). From 
Eq. (19), it is evident that systems with large flux ratios (i.e. high g) are very important 
when constraining galaxy density profiles since the flux ratio in these systems are more 
sensitive to changes in rj. If such systems are absent from the sample, the error on 77 will 
increase significantly. Also, increasing the size of the external shear by a factor of ten 
causes a systematic bias in the determination of 77 of 77 = 2 ^ 1.65. It should be noted 
that the use of flux ratio measurements is potentially very problematic because of large 
uncertainties in the microlensing magniflcation of individual images. Since the errors 
are probably not gaussian, additional information on the probability distribution of the 
microlensing magniflcations is most likely needed in order to be able to use the flux 
ratio information in an unbiased way. On the other hand, observations of the microlens 
magniflcation distribution is potentially a powerful probe of the small scale structure in 
galaxy halos. 

A possible complication is that our results for individual multiple image systems 
are based on the assumption that we can deflne average properties of lensing halos, 
i.e. the value of 7/, while N-body simulations indicate large dispersions in the individual 
halo properties [24-28]. Another caveat is that real lensing systems are too comphcated 
to allow for an analysis as simplistic as in this paper. However, we are able to constrain 



Strong lensing, cosmology and lensing halos 



13 



77 and h even with a quality factor / as low as 0.05 jj. If it is possible to retain a large 
number of lens systems, properties like redshift evolution of density profiles, variations 
in density profiles with galaxy luminosity, colour and many other important parameters 
can be measured by dividing the lens systems into different categories (in which also 
the scatter in individual halo properties may be smaller). 

Undoubtedly, modelling of real future lens systems will use more complicated lens 
models and observational information than the one used in this paper. Nevertheless, 
we believe that our simple model is able to show the general parameter dependencies 
and approximate confidence contours that can be obtained with future strong lensing 
data. Since our results are fairly robust to changes in the quantity and quality 
of the observational data, we conclude that strong lensing should be able to give 
useful constraints on the matter distribution in galaxies and galaxy clusters as well as 
complement other cosmological probes of the Hubble parameter and the matter density. 
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Appendix A. Ellipticity and external shear 

We modify our simple simulated lensing systems by adding the effects of ellipticities and 
shear. In order to treat possible non-sphericity of the lens systems we add a quadrupole 
term to the projected potential, ^, 



where is the spherical potential given in Eq. (7) and is the angle of the image relative 
to the quadrupole axis. If the system is assumed to have relatively small ellipticity, the 
angle of the second image can be written as a function of the first 



where 5 is a small parameter (and explicitly zero for spherical systems). 

Most lens systems arc embedded within an external potential which gives rise to 
an additional shear contribution. We write the combined potential as 



where — 0o is the angle between the image position and the quadrupole axis of the 
external potential. Since 7 should be small in order to justify treating the external 

potential as a quadrupole, wc derive expressions for At and r which arc valid to second 
order in q, 6 and 7 for arbitrary values of 0o- The expressions are given in Ref. [52]. 
The second order Taylor expanded expressions are excellent approximations for small q 
and/or 77 ~ 2. In practice however, we have used the full expressions for the time delays 
and the flux ratios. Since 77 ~ 2 in our simulations, the use of Taylor expansions gives 
close to identical results. 




(A.l) 



02 = 01 + TT - 5, 



(A.2) 



* = *o (1 + a cos 20) + 7extr^ cos[2(0 - 0o)], 



(A.3) 



